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Abstract. Explanations to the accelerated expansion of the Universe are usually 
sought cither in modifications of Einstein gravity or in new forms of energy density. 
An example of modified gravity is the braneworld Dvali-Gabadadze-Porrati (DGP) 
model which is characterized by a length scale which marks the cross-over between 
physics occurring in our four-dimensional brane and in a five-dimensional bulk space. 
An example of dark energy is Chaplygin gas which has similar asymptotic properties 
at early and late cosmic times. Since Chaplygin gas gives too much acceleration we 
combine it with the self-decelerating branch of the DGP model, taking the cross- 
over scales to be proportional. This 3-parameter model fits supernovae data with a 
goodness-of-fit equalling that of the ACDM model. In contrast to generalized DGP 
models and Chaplygin gas models, this model is unique in the sense that it does not 
reduce to ACDM for any choice of parameters. 
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1. Introduction 

The observation that the Universe is undergoing an accelerated expansion has stimulated 
a vigorous search for models to explain this unexpected fact. In General Relativity the 
Einstein tensor G^i, encodes the geometry of the Universe, the stress-energy tensor 
T^jy the energy density. Thus modifications to G^u imply alternative geometries, 
modifications to T^^ involve new forms of energy densities, called dark energy. The 
traditional solution is the cosmological constant A which can be interpreted either as a 
modification of the geometry or as a vacuum energy term in Tn„. This fits observational 
data well, in fact no competing model does better. 

A simple and well-studied model of modified gravity is the Dvali-Gabadadze- 
Porrati (DGP) braneworld model (Dvali & al. [T], Deffayet & al. [2]) in which our four- 
dimensional world is an FRW brane embedded in a five-dimensional Minkowski bulk. 
The model is characterized by a cross-over length scale Vc such that gravity is a four- 
dimensional theory at scales r -C where matter behaves as pressureless dust. In the 
self-accelerating branch, gravity "leaks out" into the bulk at scales r ^ Tc and the 
model approaches the behavior of a cosmological constant. To explain the accelerated 
expansion which is of recent date {z ~ 0.5), Tc must be of the order of Hq^. In the 
self-decelerating DGP branch, gravity "leaks in" from the bulk at scales r ^ Vc- Note 
that the self-accelerating branch has a ghost, whereas the self-decelerating branch is 
ghost-free. 

A simple and well-studied model of dark energy introduces into T^,^ the density 
and pressure of an ideal fluid called Chaplygin gas [31 H] . Like the previous model 
it is also characterized by a cross-over length scale below which the gas behaves as 
pressureless dust, at late times approaching the behavior of a cosmological constant. 

Since the standard 2-parametric Chaplygin gas model causes too much acceleration, 
we propose to combine it with the standard 2-parametric self- decelerating DGP model, 
taking the cross-over length scales to be proportional. The proportionality constant 
subsequently disappears because of a normalizing condition at z = 0. Thus the model 
has only one parameter more than the standard ACDM model. 

Generalizations of the Chaplygin gas model and the DGP model also have at least 
one parameter more than ACDM, yet they fit data best in the limit where they reduce 
to ACDM. This model is a unique alternative in the sense that it does not reduce to the 
ACDM model for any choice of parameters. 

2. Cross-over length scales 

On the four-dimensional brane in the DGP model, the action of gravity is proportional to 
Mpj whereas in the bulk it is proportional to the corresponding quantity in 5 dimensions, 
M|. The cross-over length scale is defined [2] as = M|j/2M|, and an effective density 
parameter as = (2rc-ffo)~^- The Friedmann-Lemaitre equation in the DGP model 
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may be written [2] 



k 1 

where a = (1 + z)~^, k = SnG/S, and p is the total cosmic fluid energy density p 
with components pm for baryonic and dark matter, and p^ for whatever additional dark 
energy may be present, in our case the Chaplygin gas. Clearly the standard FLRW 
cosmology is recovered in the limit rc ^ oo. In the following we shall only consider 
k = flat-space geometry. The self-accelerating branch corresponds to e = +1, we shall 
consider the self-decelerating branch with e = — 1. The free parameters in the DGP 
model are then and Qm = i^Pm/Ho- 

The Chaplygin gas has the barotropic equation of state = —A/p^ [31 H], where 
y4 is a constant with the dimensions of energy density squared. The continuity equation, 
p^ = —3H{p^ — A/p^), integrates to 

p^{a) = ^/A + B/a^ , (2) 
where B is an integration constant. The limiting behavior of the energy density is 

p^{a) oc VB/a^ for a ^ p^{a) oc for a > {B/AY^^. (3) 

We choose the two scales, and (B/AY^^, to be proportional by a factor H^, so that 

fV".,^;.^. (4) 

The notation is chosen because this parameter has the dimension of Hq. Note 
however, that Hq = H(l) is a constant in Eq. ([T]) (as well as in the standard model), 
whereas Hq is a function of the parameters entering the combined model. 

It is now convenient to replace the parameters A and B in Eq. ^ by Qa = k^/A/ Hq 
and H^ = Hq^W^XB / Af/^ . The dark energy density can then be written 



p^{a) = HIk-^^a^I + (i7o7^o)^(4a..«') • (5) 
3. The combined model 

In the DGP model the expansion history if (a) is given by Eq. ([T]). To distinguish our 
combined model (in which e = —1), we denote the expansion history if^(a). Making 
use of the deflnitions of Vtr^, = ^m/o^, ^a, Hq, and p^{a), the Friedmann equation 
becomes 

ii^(a)_ r„ , „ - 



+ nla-' + ^A^/l + {H^Q/HQnAn,a^) -3 .(6) 

Hq L J 

Note that VL^^ and VLa do not evolve with a, just like VL\ in the the ACDM model. 

A closer inspection of Eq. ([6]) reveals that it is not properly normalized at a = 1 
to a constant H'^{1)/Hq, because the right-hand-side takes different values at different 
points in the space of the parameters VL^, VLj-^^VLa, and Hq. This gives us a condition: 
at a = 1 we require that H'^il) = Hq. 
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This condition is a 6:th order algebraic equation in the variable Hq/Hq. Finding 
the real, positive roots and inserting the smallest one in Eq. ([6]) normalizes the 
equation properly, and rids us of the parameter Hq, which henceforth is a function 
Hq = /(f^m) ^r^,^A)- The only problem is that the function cannot be expressed in 
closed form, so one has to resort to numerical iterations. 

In this brief presentation we shall make the approximation H^{1) = Hq which is 
very good, except for at z <^ 0.1, leaving the general solution to a later paper (in 
preparation). The correct value of Hq/Hq is about 0.96. 

4. Data, analysis, and fits 

The data we use to test this model are the same 192 SNela data as in the compilation 
used by Davis & al. [3] which is a combination of the "passed" set in Table 9 of Wood- 
Vasey & al. [6J and the " Gold" set in Table 6 of Riess & al. [7] . 

We are sceptical about using CMB and BAO power spectra, because they have 
been derived in FRW geometry, not in five-dimensional brane geometry. The SNela 
data are, however, robust in our analysis, since the distance moduli are derived from 
light curve shapes and fluxes, that do not depend on the choice of cosmological models. 

We compute luminosity distances and magnitudes for all the supernovae using the 
expression ([6|. The sum turns out to be very insensitive to values of and far 
from the best fit, a problem we cure by adding a weak CMB prior, = 0.239 ± 0.09. 

The calculations are done with the classical CERN program MINUIT (James 
& Roos[8j) which delivers xiest^ parameter errors, error contours and parameter 
correlations. We do not marginalize, but quote the full, simultaneous confidence region: 
a la error contour corresponds to xlest + 3.54. The best fit is at the parameter values 

nl^ = 0.26 ± 0.16, = 0.82+°:^^, = 2.21+°:^° (7) 

with = 195.5 for 193 - 3 d.f. {x^/d.f. = 1.029), exactly the goodness-of-fit of the 
AC DM model for the same data. 

In Fig. 1 we plot the best fit confidence region in the vs Q^-plane, a banana- 
shaped closed contour. A cross in the Figure marks the point of best fit ([t]). 

A further constraint is obtained from U/r chronometry by Frebel & al. |9] of the 
age of the oldest star HE 1523-0901, t = 13.4 ± 0.8 (stat) ± 1.8 (U production ratio). 
This can be written as a lower limit to the age of the Universe t{Universe) > 12 Gyr 
(68%C.L.). The resulting forbidden region is blue in Fig. 1. The best point obtained 
from the supernova fit is now in the blue region, thus this age condition has a real 
infiuence on the allowed parameter space. 

5. Conclusions 

We embed standard Chaplygin gas in the ghost-free self-decelerating DGP geometry 
with an extra condition: that the cross-over length scales, and (B/Ay^^ respectively. 
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are proportional. The proportionality factor is fixed by a normalizing condition at 2; = 0, 
leaving three free parameters to be adjusted - one more than in the standard ACDM 
model. The resulting cosmic acceleration fits supernova data [5], El [7] excellently. It 
can also be constrained by the distance to the Last Scattering Surface, and by a stellar 
lower limit to the age of the Universe |H]. 

Generalizations of the Chaplygin gas model and the DGP model also have at least 
one parameter more than ACDM, yet they fit data best in the limit where they reduce 
to ACDM. This model is a unique alternative in the sense that it does not reduce to the 
ACDM model for any choice of parameters. 

The effective equation-of-state changed from super-acceleration to acceleration 
sometime in the range < z < 1, approaching Weff = — 1 in the future. The 'coincidence 
problem' is a consequence of the time-independent value of Tc, a braneworld property. 




Figure 1. The closed contour is the confidence region in the (fij^, rir.^)-plane from 
a fit to SNela data. The point of best fit is marked by a cross. The blue region is 
forbidden because it would make th Universe too young. 
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